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Abstract 

This paper discusses the empirical validity of Goodwin's (1967) macroeconomic model of growth with cycles by 
assuming that the individual income distribution of the Brazilian society is described by the Gompertz-Pareto distri- 
bution (GPD). This is formed by the combination of the Gompertz curve, representing the overwhelming majority of 
the population (~ 99%), with the Pareto power law, representing the tiny richest part (~ 1%). In line with Goodwin's 
original model, we identify the Gompertzian part with the workers and the Paretian component with the class of cap- 
italists. Since the GPD parameters are obtained for each year and the Goodwin macroeconomics is a time evolving 
model, we use previously determined, and further extended here, Brazilian GPD parameters, as well as unemployment 
data, to study the time evolution of these quantities in Brazil from 1981 to 2009 by means of the Goodwin dynamics. 
This is done in the original Goodwin model and an extension advanced by Desai et al. (2006). As far as Brazilian 
data is concerned, our results show partial qualitative and quantitative agreement with both models in the studied 
time period, although the original one provides better data fit. Nevertheless, both models fall short of a good empir- 
ical agreement as they predict single center cycles which were not found in the data. We discuss the specific points 
where the Goodwin dynamics must be improved in order to provide a more realistic representation of the dynamics 
of economic systems. 

Keywords: Income distribution; Pareto power law; Gompertz curve; Brazil's income data; Goodwin model; 
Growth-cycle macroeconomics; Fractals 



1. Introduction 

It has been noted long ago by Karl Marx that capitalist production grows on cycles of booms and busts. During a 
boom, profits increase and unemployment decreases since the workers are able to get better jobs and higher salaries 
due to shortage of manpower to feed the growing production. However, this boom is followed by a bust since less 
unemployment reduces the profit margin, whose recovery is achieved by a higher unemployment and a reduction of 
the workers' bargaining power. Smaller salaries increase the profit margin leading to renewed investment and then a 
new boom starts, being followed by another bust, and so on 

U Chap. 25, Sect. 1]. 

A century later Richard Goodwin [2] proposed a mathematical model which attempts to capture the essence of 
Marx's dynamics described above. In this model the basic dynamics of a capitalist society, as qualitatively described 
by Marx, is modeled by means of a modified Lotka- Volterra model where predator and prey are represented by workers 
and capitalists. Goodwin replaced the classic Lotka- Volterra dynamics of number of predators and preys by two new 
variables u and v, the former giving the workers' share of total production, which is an indirect way of describing 
the profit margin of capitalists, and v representing the employment rate, which is an indirect way of describing the 
share of those marginalized by the production, the unemployed workers, that is, the industrial reserve army of labor 
in Marx's terminology. In a boom the employment rate v increases and the workers' share u starts to increase after 
a time lag, meaning a decrease in profit margin. When employment rate is at its maximum this corresponds to the 
lowest profit margin, then the burst phase starts with a decrease in v. At this point u had already started diminishing. 
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The essence of the model is captured as a closed orbit in the u-v phase space. Clearly these two variables are out of 
phase in time J3, PP- 458-464]. 

Although the brief description given above appears to indicate that Goodwin was able to capture Marx's obser- 
vations, the model has in fact several shortcomings, the most severe one being its inability to predict quantitatively 
the above described dynamics (see below). The model was presented simply as an heuristic reasoning capable of 
giving a mathematical dressing to Marx's ideas. It was born out as a vision of the world rather than from a real-world 
data-inspired model in a physical sense. Despite this, or, perhaps, because of this, since its formulation Goodwin's 
model has attracted considerable theoretical attention in some economic circles and several variations of the original 
model were proposed [see|4|4l8l and references therein]. 

However, interestingly enough, almost half a century after its proposal, attempts to actually test this model em- 
pirically are still extremely limited. Although Goodwin's growth-cycle model is certainly influential in view of the 
number of theoretical follow-up papers cited above, studies seeking to establish its empirical soundness are limited 
only to Refs. This is a surprisingly short list when we consider the time span since the model's initial 

proposal. So little interest in empirically checking the model, especially among those who appear to have been se- 
duced by its conceptual aspects, is even more surprising if we bear in mind that for the last 30 years or so we have 
been living in an era where large economic databases are easily available digitally, so large-scale checking of this 
model against empirical data ceased long ago to pose an insurmountable barrier. Besides, even the very few studies 
which actually attempted that, all point to severe empirical limitations of the model, ranging from partial qualitative 
acceptance to total quantitative rejection. From an econophysics viewpoint it is curious that a model with such a poor 
empirical record became so influential. 

Despite this, the model does have some general empirical correspondence to reality on a qualitative level and 
this justifies further empirical studies with different databases, data handling methods and/or data type approaches. 
The basic aim must lie in identifying as clearly as possible where the model performs poorly in order to propose 
amendments and modifications. Any model, especially those theoretically seducing, can only remain of interest if it 
passes the test of experience, if it survives confronting its predictions with empirical data. If it does not survive this 
test the model must be modified, or abandoned. 

This paper seeks to perform an empirical study of the Goodwin growth-cycle model using individual income data 
of Brazil. The study presented here was directly motivated by our previous experience in modeling Brazil's income 
distribution, whose results suggested a Goodwin type oscillation in the share of the two income classes detected in 
the data 1 261 271. Building upon our previous experience with this database, we obtained yearly values of the two 
main variables of the Goodwin model, the labor share u and the employment rate v. Nevertheless, differently from 
all previous approaches for testing Goodwin's model, here the labor share was obtained by modeling the individual 
income distribution data with the Gompertz-Pareto distribution (GPD) and identifying u with the Gompertzian, less 
wealthy, part of the distribution fl27ll . The employment rate was also estimated from the same database, that is, from 
Brazil's income distribution, using the concept of effective unemployment. 

We show that from 1981 to 2009 u and v do cycle in a form bearing similarities to what the Goodwin model 
predicts, that is, closed cycles. However, our results show the absence of a single cycling center and also are in 
complete disagreement with the ones for Brazil as reported by Ref. 1 25], whose analysis employed Harvie's method 
i22ll . In addition, we attempted to see if our findings bring empirical support to the Desai-Henry-Mosley-Pemberton 
(DHMP) extension of the original model |@]. Our results show that this particular variation of the Goodwin dynamics 
has some empirical soundness, although it provides a somewhat poorer data fit as compared to the original model and 
also leaves three parameters to be determined by other, still unknown, means than the ones studied here, whereas the 
original model leaves two parameters in a similar situation. We conclude that these two models provide partial qual- 
itative and quantitative agreement with real data, at least as far as empirical data from Brazil are concerned, but both 
of them, and perhaps all variations of the original Goodwin growth-cycle dynamics, require important modifications 
and amendments before they can be considered viable representations of the real dynamics of economic systems. 

The paper is organized as follows. Sect. |2]presents a brief review of the original Goodwin model and its DHMP 
extension, focusing mostly on their dynamical equations, although some discussion about the underlying economic 
hypotheses and foundations of the original model is also presented. In Sect. [3] after a short discussion about methodol- 
ogy, we review the main equations behind the GPD. Sect. |4]analyzes the individual income data of Brazil and presents 
the u-v orbits in the 1981-2009 time period. Sect. |5]provides time variations of the employment rate as compared to 
workers' share so that line fittings allow us to determine some of the unknown parameters of both models. Finally, 
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Sect.[6]discusses the results and presents our conclusions. 



2. The Goodwin growth-cycle macro-economic dynamics 



2.1. The original growth-cycle model 

The model proposed by Goodwin is essentially a Lotka-Volterra predator-prey system of first order ordinary 
differential equations which can be written as follows 1I2II9I. I22I1 . 



ii — [— (a + d) + hv] u, 



100 - u 



(a + b) 



(1) 



(2) 



where the dot denotes the time differentiation d/df. The five constants a, b, c, d, h come from the economic hypotheses 
of the model and are supposed to obey the following conditions 1I3 H22I1 . 



c> 0, 

h >0, 

(a + d) > 0, 

(a + b)c< 100. 



(3) 



The solution of equations ([TJ and (|2j produces a family of closed orbits with period T, all having the point (u c , v c ) as 
their unique center, according to the following equations 112211 . 



: 100 - (fl + b)c, 
(a + d)/h, 



(4) 



T = 2tt I V(a + d)[l00/c - (a + &)]. 



Since u is the percentage share of labor, or workers, in national income and v represents the proportion of labor force 
employed, they both should lie in the [0, 100] interval. Here we follow the normalization adopted in Refs. 11261.12711 and 
shall refer to the maximum share, or proportion, as 100%. The upper singular point v s for the employment proportion 
is reached when v = 0, then u s = 100 - c(a + b). Similarly, when u — we have v s = (a + d)/h. However, if (a + b) is 
negative, then u s > 100, which, in principle, should not be allowed (for a conceptually possible, but so far untested, 
exception, see Ref. J2I], p. 461). Similarly, it is possible to have v s > 100. 

In this model u represents the population density of predators whereas v represents the prey population density. 
This can be seen as follows. When u = 0, ii — and v > 0. In other words, u remains equal to zero whereas v grows 
without bound, a situation happening to the prey population v in the absence of predators u. On the other hand, when 
v = 0, equations ([TJ and © together with conditions (0 show that v = and u < 0. So, without prey (v = 0), the 
predator population decreases (it < 0). 

The model is defined in terms of five parameters. However, once they are grouped as below, 



a\ — (a + d), 
a 2 = (100/e) - (a + b), 
bi =h, 
b 2 = 100/c, 

they allow equations ([TJ and (f2j to be rewritten in the form of the classical Lotka-Volterra equations fl, 

ii/u — —a\ + b\ v, 
v/v — «2 — i>i u, 



(5) 



(6) 
(7) 



that is, in terms of four parameters which could, in principle, be determined observationally, provided that both 
variables and their derivatives are obtained from real data. 
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2.2. The Desai-Henry-Mosley-Pemberton (DHMP) extension 

Desai et al. |@1 noted that the original Goodwin model can produce solutions outside the u-v domain [0, 100] x 
[0, 100] because, as seen above, both u s and v s can grow above 100. This is the main reason which led them to propose 
a modified version of Goodwin's original model, dubbed here as the DHMP extension. They also relaxed two other 
economic hypotheses assumed in the original model. So, in the DHMP extension all profits are not always invested 
and the Phillips curve, relating unemployment and inflation rate, is non-linear. Thus, the final equations yield, 



it = \-{a + d) +A(100- v) s ] u, 

V = {[-/lln(100 - u) - (a + b)] + Aln(u - w)} v, 
where d, b, d, h, 6, A, u are constants obeying the following constraints, 



(8) 
(9) 



(10) 



S>0, 
A>0, 

u <u< 100, 
h <{a + d), 

(-?-)> i. 

UOO-w/ 

Ref . Jgl] gives a clear meaning to the parameter u as being "the maximum share of labor that capitalists would tolerate", 
"typically" given by the last constraint equation in the set of expressions (ITOb above. Clearly this implies that li > 50%. 
One must also note that both the original Goodwin model and its DHMP extension consider that the labor share and 
profits are not given in terms of money, but in real terms. As we shall see below, this requirement does not pose a 
problem for our approach since our variables are currency independent [see (26)]. 

As seen above, the DHMP extension of Goodwin's growth cycle model is defined by seven parameters which can 
be grouped as below, 

ci\ — (a + d), 

ci 2 = /lln(100 - u) + (d + b), 
b\ = k 
b 2 = A, 

allowing us to rewrite equations (|8} and (|9) as follows, 



(11) 



ii/u — —a\ + b\ V , 
i'l v = -«2 + k>2 In H, 



where 



and the unemployment rate given by, 



V = 100 - v. 



(12) 
(13) 



(14) 



(15) 



Although the basic motivation for the DHMP extension was to avoid the variables of the model having values 
above 100%, this difficulty can be avoided if both u and v are defined by real data, in which case the desired threshold 
will be achieved by construction. Besides, the DHMP model has the additional disadvantage of requiring seven, rather 
than five, unknown parameters. 
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2.3. Interpretation of the conflicting variables 

As seen above the Goodwin model is essentially a predator - prey type one and this means that its two variables 
represent the opposing, but interdependent, nature of a predator - prey conflict. This is the reason why this model 
is also known as "Goodwin's class struggle model." The nature of this "struggle" arises from the possible ways we 
interpret its variables. 

On one hand, the employment rate v can be identified with the workers' class and the profit share of the "capital- 
ists" is then given by, 

UslOO-u. (16) 

In other words, U is the share of total national income obtained by the class that controls the capital, the investors. In 
this case the conflict is between the workers and the investors (capitalists). That can be seen in the light of a change of 
variables such that when u - 0, u - 0, U = and v > 0, meaning that when the profit U attained by investors remains 
constant, i.e., U — 0, the workers' share v grows without bound and represents the prey, whereas the investors U are 
in the role of predators. Here U is assumed to have a maximum value equal to 100%. 

On the other hand, following Solow 12111 . employed workers can be identified with the workers' share u and 
unemployed workers with the variable V. In this case the conflict is between employed and unemployed workers. 
When u — 0, u = and V < 0. This is consistent with the employed workforce u in the role of prey, the unemployed 
workers V being identified with the predators and the investors as passive non-players. 

However, these interpretations should not be taken at their face values as they are dependent on the conditions 
given by equations (0. Such parameter constraints were, however, not established from an analysis of real-world 
data, but came from entirely heuristic, and so far very poorly tested, reasoning. In addition, since as seen above the 
variables u and v can be identified in more than one way, this means that such interpretations must be done with 
care and always in the light of real-world data analysis and not on a speculative basis. As further emphasis of these 
difficulties, one may even argue that the constants of the model may not be constants at all, but time dependent 
variables themselves (see below). 

2.4. Origins of the Goodwin model 

As noted above, when developing his model Goodwin aimed at putting in mathematical form Marx's conceptual 
ideas about cycles in capitalism. However, as pointed out by Keen [28], Goodwin also wished to show how cyclical 
behavior could arise from very simple economic hypotheses. Next we shall present a simple derivation of the model 
in order to highlight that it results from an extremely simplified representation of the economy. 

Let K be the amount affixed capital (plant and equipment) and Y the output that an economy can generate. The 
output to capital ratio cr clearly varies over time in a country, but let us consider it a constant as a first approximation 
and write it as follows, 

cr = K/Y. (17) 

If L is the amount of labor for a given output, one can also assume as first approximation a constant output to labor 
ratio a, that is, 

a = Y/L. (18) 
The amount of labor can be written in terms of the population N and the employment rate v as follows, 

L = Nv. (19) 

Let w be the average wage value. Then the wage bill, that is, the total amount of wages in an economy is given by, 

W = wL. (20) 

At a first approximation the employment rate can be related to the rise of wages as follows, 

w/w = /i(v). (21) 

Since the wage share u is given by, 

u = wL/Y, (22) 
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remembering that a is constant, equation ( f2TT > becomes, 



m = fM. (23) 

This expression reduces to equation © if f\ (v) is assumed to be a linear function. 
The profit level P is given by, 

P = Y - W. (24) 
As a first approximation all profits are invested, so the profit share P/ Y is the investment Y. Hence, 

T = P/Y= 100 -(W/F) = 100 -k = U. (25) 

Here the unit was changed to 100% due to our previous choice of normalization. The profit rate n is given by, 

„ = p/K, => T = cm — 100 — u, (26) 

which can be rewritten in functional form as below, 

Y = f 2 (u). (27) 
Investment is also the rate of change of capital K/K. So, 

K Y v 

T= — = - = - + const. , (28) 
K Y v 

where the constant comes from the hypothesis of a steady labor supply, e.g., L changes exponentially. Summing up 
we have that, 

v/v = / 2 (w), (29) 

which reduces to equation © if fz(u) is assumed linear. 

Clearly the model results from extremely simple specifications of the economy. But, it is so simple that it cannot 
reproduce the frequency properties of output growth in a certain time period or the distribution of recession sizes and 
duration. However, the dynamic stochastic general equilibrium (DSGE) models of cycles adopted by current neoclas- 
sical economics cannot do so either M29H31I1 . hence what is remarkable is that the very restricted model proposed by 
Goodwin finds any empirical support in real data [28]. 



3. The Gompertz-Pareto income distribution 

Econophysics is a new research field whose problems interest both economists and physicists. However, when 
physicists approach a problem traditionally dealt with by economists, they do so under a very different modeling 
perspective. Although it is uncommon to find methodological issues discussed in physics papers, considering the 
hybrid nature of econophysics and the theoretical crisis of the current mainstream economic thought Il29l442ll . it is 
worthwhile to emphasize the differences in methodological perspectives between physics and economics regarding 
model building and, especially, model abandoning. We have already expressed some of our thoughts on this topic in 
Ref. 11261 Sect. 3], but a few more words are worth saying before we review our approach to the income distribution 
problem. 



Econophysics was born and remains a branch of physics I143l - l45ll . employing, therefore, its centuries old proven 
epistemological methodology. It considers a scientific theory as being made by laws of nature, which are theoretical 
constructs, often expressed in mathematical language, that capture regularities, processes, structures and interrelation- 
ships of reality. Successful physical laws provide good empirical representations, or images, of the real world, of 
nature, and allow us to reach predictions regarding the outcomes of processes that do go on in nature. However, by 
being images of nature, these laws are obviously limited and, hence, they will always provide imperfect representa- 
tions. The only way we can ascertain how imperfect they are is by practice, i.e., by creating pragmatic measures of 
the adequacies of these laws, always empirically comparing their predictions with what occurs in the real world [46 j. 
In other words, good laws provide good predictions, bad laws provide bad predictions. This has nothing to do with 
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the extensive use of mathematics by physical theories. Mathematics is a language, a tool of formal logic, and by itself 
has no a priori relationship with physical, or social, reality. Physicists choose if and which mathematical tools are 
required to express something observed in naturefl 

Since our understanding of the theories behind these laws changes with time, the same occurring with the measures 
of adequacies due to technological advances, we must keep measuring the adequacies of these laws by perfecting old 
measures as well as creating new ones, that is, constantly updating our theories and models through practice in order 
to find their limits of validity. The theoretical aspects behind these laws, even their metaphysical presuppositions, 
must also be perfected by shedding the inappropriate elements so that the appropriate residue remains, in a process 
very similar to Darwin's natural selection. And, if there is no appropriate residue left the theoretical construct is 
abandoned, becoming extinct J49ll . Under this viewpoint, a model is a more restricted theoretical construct, taking one 
or two elements above - regularities, processes, structures and interrelationships -, but not all of them. Nevertheless, 
a model is also subject to measures of adequacy and since they incorporate less elements than a theory, it suffers a 
more rapid process of perfection by selection as well as extinction. 

Physicists have been following this methodological approach for centuries and as a consequence they have amassed 
a large number of physical theories that were perfected by generations of physicists, who kept their appropriate kernels 
but changed their original elements in various degrees, and also to many other theories which are now superseded. 
Theoretical pluralism is tacitly accepted as an essential element for the development of physics. Real science starts 
from observation of nature, either physical or social, and any theoretical discussion must keep referring back to em- 
pirics, a factor that limits and guides any theoretical debate, leading to healthy refining, replacing or even abandoning 
of theories and models ifsoh . 

However, it seems that this methodological viewpoint regarding model checking has not been adopted by a sizable 
number of economists. Econophysicists are often perplexed to witness how often economists confuse their models 
with reality, showing a behavior which was already described as 'scientific dogmatism' 114611 . Thus, they would of- 
ten disregard startling obvious empirical facts rather than change or dismiss their inappropriate theories or models 
151 . 521. showing to a large extent an absolute devotion to theoretical economic constructs, especially an empirically 



unwarranted obsession with equilibrium, in parallel to little or no empirical interest, often keeping such a theoretical 
worship even when empirical evidence that might support the theory is absent. Worse still, even when there is evi- 
dence that directly contradicts what would be predicted to occur by applying the theories 11531 pp. 2-5]. Some would 
say this phenomenon is due to 'ideological assumptions', disguised visions of the world under scientific pretenses 
B48I1 . Others call this behavioral mode 'cargo cult economics' 15411 in reference to the famous Feynman speech about 
methodologically inadequate, or false science J55, 56]. Nevertheless, the epistemological ideas above, adopted by 



physicists a long time ago, are apparently being slowly absorbed into the economic thought 0571, 15811 . 



Having stated our methodological viewpoints, next we shall review the basic hypotheses and equations behind the 



GPD as advanced in Refs. ||26L|27|1 . 



3.1. Definitions 

Let T{x) be the cumulative income distribution giving the probability that an individual receives an income less 
than or equal to x. Then the complementary cumulative income distribution F(x) gives the probability that an individ- 
ual receives an income equal to or greater than x. It then follows that T{x) and F(x) are related as follows, 

T(x) + F(x) = 100, (30) 

where the maximum probability is taken to be 100%. Here x is a normalized income obtained by dividing the nominal 
income values by some suitable nominal income average l26ll . If both functions T{x) and F{x) are continuous and 



'Here we take a viewpoint different from Lawson's 14711 regarding the role of mathematics in economics, a viewpoint based on the larger 
experience of other sciences which successfully adopted mathematical modeling, especially, but not restricted, to physics. The obvious failures 
of mathematical modeling in economics is a problem specific to academic economics because it misinterpreted the role of theoretical thinking by 
means of a continuing excessive emphasis in theoretical introspection parallel to a strong downplaying of the empirical certification of models. 
Hudson 1 48] provides an interesting account of why and how academic economics reached this present state of affairs. One must note that the 
impressive achievements of the 20th century in theoretical physics would never had occurred if physicists had ignored empirics to the extent that 
academic economists do. 
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have continuous derivatives for all values of x, we have that, 



dT(x)/dx = f{x), dF(x)/dx = -f(x), 



(31) 



and 



f 

Jo 



f(x) : dx = 100, 



(32) 



where f(x) is the probability density function of individual income. Thus, f(x), dx is the fraction of individuals with 
income between x and x + dx. The equations above lead to the following results, 



fix) - <F(0) 
whose boundary conditions are, 



I f(w) dw, 
Jo 



F(x) - F(oo) 



f 



f(w) dw, 



T(0) = F(oo) =s 0, 

r(oo) = f(0) = ioo. 



(33) 



(34) 



Clearly both T(x) and F(x) vary from to 100. 

3.2. The Gompertz-Pareto distribution (GPD) 

The GPD was proposed in Ref.JIi] and discussed in detail in Ref. [27] . Its complementary cumulative distribution 
is formed by the combination of two functions which can be identified with the two main classes forming most modern 
societies, workers and investors (capitalists). The first component describes the lower part of the distribution, that is, 
those who survive solely on their wages, the workers, and is given by a Gompertz curve. The second component of 
the complementary cumulative distribution describes the tail of the distribution by means of the Pareto power law and 
represents the investors, that is, the rich capitalists. Then we have that, 



F(x) = 



G(x) 



JA-Bx) 



*» ^ A - B ^ x -a 



{ P(x) = (x,) a e 1 
and the cumulative income distribution may be written as below, 

JA-Bx) 

&(x) = 100 - e' 

T(x) 



( < x < x,), 

(X, < X < oo), 



(Gompertz) 
(Pareto) 



(35) 



P(x) = 100 - (x,) a e' 



,(A-Bx,) 



( < x < x,), 



(X, < X < oo). 



(36) 



Here x, is the income value threshold of the Pareto region, a is the Pareto index describing the slope of the power law 
tail, B is a third parameter characterizing the slope of the Gompertz curve and A is a number whose value is set by 
boundary conditions, as follows. Since G(x) = exp[exp(A - Bxj], the condition ( 1341 implies G(0) = 100, then we 
have that, 

A = In (In 100) = 1.5272. (37) 

(A-Bx t ) 

The term (x t ) a e is the normalization constant of the Pareto power law and comes as a consequence of condition 

(l32l . as well as the continuity of functions d35l > across the frontier between the Gompertz and Pareto regions, defined 
to be x — x,. 

The equations above allow us to find the expressions for the probability density income distribution, 

JA-Bx) 

g(x) = B e^"- 1 ' e' 

/(*) = 



p(x) = a (x,) e c x 



( < x < x t ), 

(x t < X < oo), 



(38) 
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as well as the average income of the whole population described by the GPD, 



<*> = 



100 



x f(x) dx = 



100 



ax, e (A-Bx,) 

I(x t ) + -e c 

(a - 1) 



(39) 



where, 

I(x)= wg(w)dw= w B ) e dw. (40) 
Jo Jo 

The parameters a, x, and B are all positive and they fully characterize the GPD. However, due to convergence require- 
ments [26], the expression ( |39l for the average income is only valid if a > 1. Both a and B can be determined by 
linear data fitting since equations (l35l > can be linearized. However, x, is independently found under the constraint that 
the boundary condition d37l > is satisfied to whatever degree of precision the available data allow. 

The Lorenz curve of the GPD has its X-axis given by the cumulative income distribution function T{x), whereas 
the first-moment distribution function T\(x) defines its Y-axis. Accordingly, they can be written as follows [27], 



T(x) 



-f 

Jo 



f(w) dw 



and 



i r 

Ti(x) = — wf(w)dw = 
(x) Jo 



AA-Bx) 

100 -e e 
100-(x,)"e e 

m 

(x) ' 



(A—Bx t ) 



(0 < x < x t ), 

(X t < X < oo), 



(41) 



100 + 



a{x,) a e' 



,(A-Bx,) 



(I -a) (x) 



(0 < x < x,), 



(x t < x < oo; ). 



(42) 



Thus, T\(x) varies from to 100 as well. The Lorenz curve is usually represented in a unit square, but the normaliza- 
tion (1321 ) implies that the square where the Lorenz curve is located has area equal to 10 4 . 
The Gini coefficient under the currently adopted normalization is written as, 



Gini =1-2x10" 



r 1 (x)f(x)dx. 

Jo 



Considering now equations d38l > and d42l . the Gini coefficient has the following expression in the GPD, 



Gini = 1 - 2 x 10" 



JA-Bx) P (A-Bx,) 

(A-Bx) Q e dx+100e e 



a 2 x t e 2e 



(A-Bxt) 



{x){a-\)(\-2a) 



(43) 



(44) 



As discussed in 112711 . we can define the percentage share of the Gompertzian part of an income distribution 
described by the GPD by means of equation d42l . This quantity may then be written as follows, 



Tx(x t ) = 100- 



* f -e e 



(A-Bx,) 



(a - 1) (x) 



(45) 



Hence, we identity the percentage share of the lower income strata described by the GPD with Goodwin's labor share 
u. Note that by doing so, u no longer represents the industrial reserve army of labor, but in fact the relative surplus 
population since the latter includes not only the unemployed, but also those unable to work. Such identification 
allows the description of the Goodwin variables in terms of measurable quantities connected to different income 
classes whose empirical values can be obtained, for instance, from the Lorenz curves. This connection can be made 
clearer by the inversion of equation (PlBl . 



1 = 1 

a 



JA-Bx,) 



Mx,e 



(100-w)J(x t ) 



(46) 
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Due to the high non-linearity of this expression one can only use it to determine a if the values of u, B and x, are 
known to a very high degree of accuracy. 

The equation ( f46b links the Pareto index a to parameters which are solely determined in the Gompertzian segment 
of the distribution: the cutoff value x,, the Gompertzian percentage share u and its distribution slope B. In other 
words, equation d46l > links the income distribution of the lower and upper classes forming a society, showing clearly 
their dynamical inter-dependency. If we consider that temporal changes in the income distribution do take place, 
we can no longer consider these quantities as parameters. Some of them, or perhaps all of them, ought to be time 
dependent variables (see below). 

The GPD requires a > 0. In addition, an average income is only possible if a > 1. Considering these two 
conditions in equation (|46*T > we conclude that, 



< 



JA-Bx,) 

ux,e 



(100 -m) JO,) 



< 1 and u< 100. (47) 



Remembering equation (fTol i the last condition is equivalent to U > 0, which means that an income distribution 
described by the GPD is only possible in a system where investors have a nonzero share of the total income. 

3.3. Exponential approximation 



As shown in Refs. [26, 27], the upper part of the Gompertz curve can be approximated by an exponential and this 
allows us to take this subdivision of the Gompertz curve as representing the middle class present in most societies. In 
other words, in this approach of the income distribution characterization of a society we assume that the middle class 
is just the upper echelon of the wage labor class. Thus, for Bx > A, e~ B v < 1 and x < x, we have that, 

( G(x) * 99 + er Bx , 

I @(x) * 1 - e-**, (48) 
[ g(x) * B er Bx , 

which are already normalized to obey the boundary conditions Q41 l. If the lower stratum of a society is formed 
essentially by a very large middle class, one can in principle write all equations shown in Sect. I3.2l in terms of the 
approximations (1481 1. although in such a case we can expect a certain degree of distortion in the distribution since all 
modern societies seem to have a certain percentage of very poor people, however small this percentage may be. 



4. Cycles in the income and employment data of Brazil 

Publicly available individual income distribution data of the Brazilian population have allowed Moura Jr. and 
Ribeiro 12611 to determine the GPD parameters from 1978 to 2005 after a careful handling of the data. Chami Figueira, 
Moura Jr. and Ribeiro 02711 extended this analysis to include income data for 2006 and 2007, as well as showing how 
the GPD produces results compatible with those obtained directly from the raw data, that is, without assuming the 
GPD, with error margins up to 7%. In this work we further extend these two previous analyzes to include data for 
2008 and 2009, but disregarding the results for 1978 and 1979 due to their unreliability [27]. 

Table Q] presents the three GPD parameters B, x, and a followed by the unemployment rate [V], Gini coefficient 
and the percentage share of the Gompertzian component of the distribution. B and a were obtained by linear data 
fitting whereas x, was determined such that a linear fit would produce the boundary condition ( l37l i with discrepancy 
of about 2%. Lorenz curves were generated from the raw distribution for each year allowing the calculation of the Gini 
coefficient without assuming the GPD, denoted here as [Gini] in order to distinguish it from the one obtained assuming 
the GPD in equation (|44| ). Once x, was found it became possible determine [u] directly from the raw data, that is, 
without using equation d45l >. Similarly, [V] denotes the unemployment data without any distribution assumptions, [v] 
is obtained using equation ( TT3T > and xj is the unemployment income threshold used to calculate [V] (see below). The 
time derivatives are given by the expressions, 

M = M = £(100-[V]). (49) 

at at 
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Table 1: Data for Brazil from 1981 to 2009. Values in brackets mean that they were evaluated without using the GPD parameters, that is, 
directly from the raw data. This table contains the GPD parameters B, x, and a for the individual income, unemployment income threshold xj, 
unemployment rate [V], Gini coefficient, percentage share [h] of the Gompertzian component (workers' share), employment rate [v] and time 
derivatives of the last two quantities, as given by equations i49\ . The results from 1981 to 2007 had already appeared in [26, 27] whereas those for 
2008 and 2009, as well as the ones for employment and unemployment, are new. The time derivatives [u] and [v] were calculated numerically using 
equation j5U . Since there were no income samplings in 1991, 1994 and 2000 [see|26J, some results for these years were obtained by numerical 
interpolation. 



year 




B 






a 


-V,/ 


[V] (%) 


[Gini] 


[«] (%) 


[u] (%/year) 


[v] (%) 


[v] (%/year) 


1981 


0.342 


+ 0.016 


7.533 


2.839 


+ 0.091 


0.182 


14.8 


0.574 


87.7 




85.2 




1982 


0.342 


± 0.015 


7.473 


2.677 


± 0.042 


0.174 


14.5 


0.581 


87.2 


+ 1.08 


85.5 


-0.20 


1983 


0.330 


±0.010 


6.910 


2.636 


±0.081 


0.175 


14.5 


0.584 


85.5 


-0.04 


85.5 


-1.06 


1984 


0.332 


±0.013 


7.388 


2.839 


± 0.072 


0.170 


12.4 


0.576 


87.2 


-0.17 


87.6 


-1.32 


1985 


0.329 


±0.010 


7.490 


2.656 


± 0.093 


0.154 


11.8 


0.589 


85.8 


+0.99 


88.2 


-2.22 


1986 


0.344 


±0.013 


7.112 


2.567 


± 0.065 


0.127 


7.9 


0.580 


85.2 


-0.05 


92.1 


-1.16 


1987 


0.343 


±0.016 


7.626 


2.724 


± 0.057 


0.127 


9.5 


0.592 


85.9 


-0.08 


90.5 


+2.11 


1988 


0.324 


±0.014 


8.140 


2.874 


±0.125 


0.133 


12.1 


0.609 


85.4 


+ 1,74 


87.9 


+0, 17 


1989 


0.317 


±0.010 


7.856 


2.428 


± 0.079 


0.111 


9.9 


0.628 


82.5 


-0.23 


90.1 


-2.35 


1990 


0.335 


±0.015 


8.074 


2.636 


± 0.053 


0.099 


7.4 


0.605 


85.9 


-1.98 


92.6 


+0.48 


1991 














10.8 




86.4 


-0.57 


89.2 


+3.37 


1992 


0.364 


± 0.020 


7.635 


2.636 


± 0.063 


0.162 


14.2 


0.578 


87.0 


+ 1.18 


85.8 


+0.53 


1993 


0.330 


± 0.008 


7.674 


2.567 


± 0.042 


0.137 


11.9 


0.599 


84.1 


+ 1.01 


88.1 


-2.54 


1994 














9.1 




85.0 


-0.92 


90.9 


-2.75 


1995 


0.333 


±0.012 


7.887 


2.777 


±0.106 


0.098 


6.4 


0.596 


85.9 


-0.86 


93.6 


-0.63 


1996 


0.347 


± 0.020 


8.163 


2.749 


±0.107 


0.096 


7.8 


0.598 


86.7 


-0.12 


92.2 


+0.43 


1997 


0.338 


±0.015 


7.935 


2.617 


± 0.052 


0.099 


7.2 


0.598 


86.1 


+ 1.09 


92.8 


-0.29 


1998 


0.326 


± 0.009 


7.628 


2.677 


±0.031 


0.103 


7.3 


0.597 


84.5 


+0.08 


92.7 


+0.25 


1999 


0.331 


±0.013 


7.811 


2.777 


± 0.068 


0.107 


7.7 


0.590 


86.0 


-0.53 


92.3 


+0.55 


2000 














8.4 




85.6 


+0.40 


91.6 


+0.66 


2001 


0.335 


±0.011 


7.774 


2.724 


± 0.205 


0.122 


9.0 


0.592 


85.2 


-0.41 


91.0 


-0.24 


2002 


0.339 


±0.015 


7.878 


2.500 


±0.121 


0.123 


7.9 


0.586 


86.4 


-0.08 


92.1 


+0.01 


2003 


0.333 


± 0.009 


7.374 


2.777 


± 0.057 


0.134 


9.0 


0.579 


85.4 


-0.40 


91.0 


-1.07 


2004 


0.333 


±0.017 


8.005 


3.234 


±0.133 


0.105 


5.7 


0.582 


87.2 


-0.44 


94.3 


-0.59 


2005 


0.326 


± 0.009 


7.403 


2.839 


± 0.089 


0.118 


7.9 


0.580 


86.2 


-0.26 


92.1 


+2.06 


2006 


0.323 


±0.015 


8.078 


3.749 


±0.136 


0.125 


9.9 


0.592 


87.7 


+0.27 


90.1 


-0.29 


2007 


0.334 


± 0.009 


6.934 


2.839 


±0.104 


0.125 


7.3 


0.572 


85.7 


+0.28 


92.7 


-1.03 


2008 


0.366 


±0.011 


6.848 


2.567 


±0.051 


0.141 


7.8 


0.543 


87.2 


-0.36 


92.2 


+0.26 


2009 


0.363 


±0.010 


6.500 


2.656 


± 0.065 


0.148 


7.8 


0.539 


86.4 




92.2 
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One should note that the focus of this paper is not to discuss the adequacy of the GPD description of income 
data by comparing results obtained by assuming or not the GPD, that is, comparing Gini to [Gini] or u to [u], as this 
task was already successfully accomplished in Ref. 12711 . Our focus here is to use the GPD as a tool to partition the 
income distribution in the Gompertzian and Paretian components, identify the former with one of the variables of 
the Goodwin model and to discuss the possible dynamical implications of such a division, that is, linking the GPD 
parameters to the Goodwin model. The unemployment data appearing in Table [T]require, however, some explanation 
about how they were determined. 

Two basic facts prevented us from using official Brazilian joblessness statistics in the analysis studied here. First, 
unemployment data collection methodology changed quite substantially during the time period of this study (1981 
to 2009) and, secondly, its sampling method differs from the one used to survey income. Taken together, these 
two facts imply the lack of sample homogeneity in the whole period of this analysis, which renders it impossible 
to derive measurable quantities without introducing substantial statistical biases. Without sample homogeneity we 
cannot compare unemployment data from early and late years in the studied time interval. These difficulties can be 
avoided if unemployment is directly estimated from the income database by means of a criterion applicable to the 
entire time period of this study. The reasoning we followed to do that is described below. 

Every society produces useful energy and materials to be consumed by the people who participate in their pro- 
duction. This means that a person active in this production receives a share of those materials and energy, that is, 
a share of the total value produced by the society in a certain period of time. Income is, therefore, a flow of value 
(energy and materials) a person receives in a certain time period. Under this viewpoint, even food is part of this share. 
The unemployed is the individual who does not participate in the production and, therefore, does not receive value. 
Nevertheless, nobody can survive too long without food or a minimum amount of energy and, thus, if the individual 
survives this means that somehow this individual still has a value inflow. Such a minimum supporting value is usually 
provided by family or, in more limited ways, by the state, but actually means a reduced value inflow for the group 
family this individual belongs. In other words, when somebody becomes unemployed those close to this individual 
are the ones who suffer most because the whole family has a smaller share of value or, which is stating the same, the 
group family income decreases. So, there should be a limit in income distribution where unemployment, or underem- 
ployment, can be effectively detectable. We call this limit effective unemployment. An average person who receives 
up to this minimum income barely participates in the production and for all practical effects is jobless. 

Following this reasoning we then probed the data for income values which would produce unemployment rates in 
agreement to those in the official unemployment surveys for the last 15 years or so. Our results showed that effective 
unemployment occurs when the average individual income is equal to or below 20% of the national minimum salary 
in Brazil expressed in US dollars and effective at the time the income survey was carried out (September of each year). 
This amount defines the unemployment income threshold xj which, after being normalized to become a currency free 
quantity, was applied to the income distribution of each year to obtain the percentage share of those in the distribution 
whose income were equal to or below this amount. This method provides our effective definition of unemployment. 

Connecting the unemployment income threshold Xd to minimum salary has the advantage of providing a simple 
criterion applicable to income data for all years of this study, even before 1994 when Brazil sampled unemployment 
through a different methodology and experienced runaway inflation and hyperinflation. The results for the unemploy- 
ment rate [V] obtained using this criterion are presented in Table Q] Note that once xj is known, the GPD allows us 
to obtain V by means of an expression similar to equation (l45l l. Indeed, as one should have Xd < x t , remembering 
equation (|42T > we conclude that, 

I{x d ) 

V = Ti{x d ) = ±^-. (50) 
(x) 

We can now plot the results for [u] and [v]. Figure Q] shows the time evolution of these two variables where one 
can see that both variables cycle with similar periods of about 4 years. In addition, these cycles are apparently out-of- 
phase for most of the studied time interval and have phase difference of about 2 years. This clearly implies short term 
cycles. These results bring qualitative empirical support to the Goodwin approach for describing the dynamics of the 
capitalist production as described by Marx. 

Figure [2] shows the [m]-[v] phase space where one can see clockwise orbits for most of the time interval, a fact 
which again brings qualitative empirical support to the Goodwin model at least as far as Brazilian data is concerned. 
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However, the orbits clearly do not have a single center, as both the original Goodwin model and its DHMP extension 
predict. After 1994 the center of the orbits seems to move to an upper position in the phase space. In order to better 
appreciate this change, figure|3]shows the same results of figure[2] but divided in two time intervals, from 1981 to 1994 
and 1995 to 2009. These results clearly contradict the Goodwin prediction of all orbital centers having the same fixed 
coordinates u c and v c , as described by equation ©. One should also note that these results are entirely different from 
the ones obtained for Brazil in Ref. [25] using Harvie's method (22|| and in much better agreement at a qualitative level 
with the Goodwin model. Finally, figure [4] shows the same data in a 3-dimensional plot with the Z-axis representing 
the time. This graph provides a different way of seeing the displacement of the points to a different region after 1994 
by means of their projection in the YZ-plane, as well as a possible earlier displacement, whose transition occurred 
from 1981 to 1983. 




1980 1984 1988 1992 1996 2000 2004 2008 

year 



Figure 1 : Time evolution of the Gompertzian component (workers' share) [«] and employment rate [v] in Brazil. The plot shows that these variables 
cycle out-of-phase for most of the studied time interval with periods of about 4 years in both variables, meaning that booms and busts in Brazil 
occur in short term cycles. These results show that predator-prey like models can be used to represent real economic systems. 



The important event which may explain why the apparent orbital center changes location after 1995 is the end of 
hyperinflation. In 1994 Brazil established a new and stable currency, the real (R$), which abruptly ended the strong 
inflationary period of the previous 15 years. This fact seems to be reflected in the [m]-[v] phase space by a change in 
the center of the orbit. One can also see in tableQ~]a slow, although modest, decrease in the Gini coefficient after 1993. 
In addition, since the Brazilian high inflationary period started at about 1980, the positions corresponding to the years 
of 1981 to 1983 in the phase space appear to represent an earlier transition from yet another region in the phase space. 
This seems to be the case if we carefully look at these points in the graphs of figures|2]and|4] 

The absence of a single center for all orbits means that the parameters a, b, c, d, and h of the Goodwin model are 
most likely not constants at all, but time dependent variables. Nevertheless, at a qualitative level the model certainly 
has empirical support which justifies the identification of [u] and [v] with u and v, although in order to understand the 
real dynamics behind these quantities one probably needs to somehow modify the dynamical equations (Q3 and (0 to 
reflect these empirical evidences. 

Finally, we should note that the lack of a single orbital center in real-world data has already appeared in earlier 
empirical studies earned out by other authors on the Goodwin model. The u-v phase-space plots of Desai 12011 . Solow 
12111 . Harvie I22I1 . More no 1I23I1 . Mohun & Veneziani [24] and Garcia Molina & Herrera Medina ll25ll show similar 
results as ours. Vadasz [10] also reached a similar conclusion, although by indirect means. The important point is 
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Figure 2: [tt]-[v] phase space for Brazil from 1981 to 2009. The plot points are labeled in growing numerical sequence with each integer number 
representing one year in the studied time interval. Thus, the label given by the number '1' indicates the year 1981 and label '29' means the year 
2009, providing then a clear visualization of the largely clockwise evolution of the cycles. One can see that before 1995, indicated by label '15', 
the system was cycling in a different region of the phase space. The end of hyperinflation in 1994 (label '14') is possibly the event which made the 
system move to a new cycling region, where it still remains. 



that all these authors reached the same conclusion despite their use of very different methods to analyze observational 
data. Therefore, one feels justified to conclude that this feature appears to be universal and clearly indicate that the 
Goodwin model must be changed in order to accommodate this real-world feature. 

5. Temporal variation of the employment rate and workers' share 

The data presented in tableQ]allow us to go beyond the qualitative discussion of the previous section and carry out 
a quantitative evaluation of the Goodwin model and its DHMP extension. To do so we need first to carry out simple 
numerical estimations of the time derivative [u]. This task is most straightforwardly accomplished using the following 
expression, 

. [«](* + AO- [u](t- At) 

* TKt ' (51) 

where At = 1 year. Similar procedure is used to determine [v]. The goal here is to use data fitting to estimate the 
parameters of the two sets of dynamical equations, the first set being given by equations © and (0 of the original 
Goodwin model and the second one by equations (IT~2b and dT3l l which constitute the DHMP extension. 

5.1. Goodwin model 

Figure [5] shows two plots, the left one for the variables of equation © and the right plot for equation (O. The 
fitted straight lines parameter values are also presented for both plots. It is clear that both sets of points are compatible 
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Figure 3: These two graphs present the same data as in Fig. [2] but divided into two sets of points. The plot in the left shows the [h]-[v] phase space 
from 1981 until 1994, whereas the right plot presents the data points from 1995 to 2009. This shows even more clearly that the Brazilian economic 
system moved from one region to another in the phase space during the time interval studied here. One can also note in the left plot that the system 
was possibly moving from yet another region in the period from 1981 to 1983, since the labels '1' to '3' indicating these years appeal' to be part 
of a transition from a different region than the one where the system remained until 1994. This possible interpretation has some empirical support 
because the high inflationary period in Brazil started in about 1980. 



with a linear approximation similar to the original dynamical equations, but the parameters behave in exactly opposite 
manner from what the model predicts. While the slope of the lines predicted by equations (JSJ and O are, respectively, 
positive and negative, the results coming from Brazilian real-world data are the other way round. This is clear in 
both graphs. This result can also be seen if we use the fitted parameters to obtain conditions which the supposedly 
"constants" of the Goodwin model should obey. Doing so we conclude that the Brazilian economic dynamics studied 
here gives, 

c<0, 
h < 0, 

' (a + d)<0, { ' 

(a + b)c> 100. 

These results completely upset the parameter conditions given by equations (0, which were thought to be valid. The 
fitting also leaves two parameters yet to be determined by some yet unknown equation relating them since, as seen 
above, the orbital center and period equations (H) are clearly invalid in the Brazilian income dynamics. 

The calculated uncertainties in the fitted parameters do not change this situation, a fact which forces us to conclude 
that the economic hypotheses advanced by Goodwin to derive his model are either not applicable, partially or com- 
pletely, to the economic system studied here or they are flawed. Whatever conclusion one may choose, this analysis 
indicates that to advance this model with the aim of turning it into a viable representation of the real world, the focus 
must lie on the probable modification of the set of differential equations (0 and (0 and their empirical validation, 
rather than how they were obtained. Only after a good model is achieved, and by good we mean a model with solid 
empirical foundations, may we start looking for the real economic conditions behind its dynamics. 

Since the data show that the parameters of the model follow the exact opposite predictions given by the expressions 
0, another consequence of the results shown in figure |5]is the reversal of the predicted roles of predator and prey 
discussed in Sect. 12.11 Indeed, according to the fitted parameters (see caption of figure [5), when [u] = 0, [u] = 
and [v] < 0. In this case [u] plays the role of prey because without it ([u] = 0) the predator population decreases 
([v] < 0). Similarly, when v = 0, [v] = and [u] > 0. So, v plays the role of predator because without them (v = 0) 
the prey population grows without bounds ([u] > 0). Such a reversal of roles of predators and preys coming from 
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Figure 4: This plot is a 3-dimensional representation of the same points appearing in figures [2]and[3] It provides a different visualization of the 
system displacements during its evolution from 1981 to 2009, showing more clearly in the YZ-plane projection (the surface plane for [v] vs. time, 
on the left side of the plot) the three regions where the system located itself in the studied time interval. The points for 1981 to 1983 (labels T to 
'3') seem to be a transition from an unspecified earlier region where the system stayed before the high inflationary period started at about 1980. 
The end of hyperinflation in 1994 (label '14') moved the system to yet another region on the top right of the YZ-plane. 



the real-world data analysis presented here also implies a reversal of the reasoning presented in Sect. [23] regarding 
how one interprets the conflicting variables. However, although such role discussions had some importance in the 
past, such interpretations are now of lesser importance than revealing the inner dynamics of these two inter-dependent 
variables. When such dynamics is better understood by means of realistic, not introspective, models, such roles will 
naturally emerge from those real-world representations. 

5.2. DHMP extension 

The variables in the dynamical equations (TT2l and (TTJb of the DHMP extension are plotted in figure [6] To do 
so we had to choose a value for the maximum share of labor [h]. From table [T] we see that the highest value in the 
studied time period is 87.7% in 1981 and in view of the fact that the DHMP model does not give any hint about how 
to obtain [u], assuming whatever constant value higher than that is enough for our purposes here and will not change 
the general behavior of equation (fT3l . So, we chose [u] = 95% as a reasonable value for this analysis. 

The left plot in figure [6] shows the points related to the dynamical variables of equation (fl2l) while the right graph 
in concerned with the variables appearing in equation (fl~3). The fitted parameters are written in the figure caption and, 
similarly to our reasoning above, they produce real-world conditions for the "constants" of the DHMP model. They 
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Figure 5: Left: [u]/[«] vs. [v]. Right: [v]/[v] vs. [it]. Although both graphs show some dispersion of the results, one can clearly identify a general 
tendency for the observational points to decrease in the left plot and to increase in the right one. Straight line fits on both sets of results, indicated as 
full lines, produced the following results. For the left plot, the expression [u]/[w] = A\ +Bi[v] resulted in Ai = 0.17 ±0.06, B\ = -0.0019 ±0.0006. 
For the right plot, the equation [v]/[v] = Ai + Bi\_u\ yielded parameters as follows: Ai = -0.52 ± 0.22, Bi = 0.006 ± 0.003. These results should 
be compared with equations |6} and 0. 



may be written as follows, 

S>0 (?), 
h > (?), 

■ (a + d)>0 (?), (53) 
A<0, 

(a + b)< -^ln(lOO-S). 

The results with a question mark are inconclusive due to the uncertainties of the fitted parameters. For the other 
two, A < contradicts the prediction given in equations (flOt . but implies that (a + b) > for the chosen u. So, 
despite the fitting, the DHMP model remains in a very inconclusive status regarding the empirical behavior of its 
dynamical variables and its supposedly constant parameters. Even so, because the model has too many parameters, 
after a successful fitting where one of the parameters had to be assumed (u), two other parameters remained unknown 
and still require determination by at least another, also unknown, expression. 

In conclusion, because the DHMP extension has more unknown quantities and its dynamics is described by some- 
what more complex differential equations than the original Goodwin model, comparing its predictions with the Brazil- 
ian data renders mixed and inconclusive results. Adding to this situation are the high errors in the fitted parameters 
and the fact that even after a successful fit several parameters remain unknown. These results place the DHMP exten- 
sion in a much less favorable situation than the original Goodwin model regarding empirical validity, at least as far as 
Brazilian data is concerned. 



6. Conclusions 

In this paper we have studied the empirical validity of the model of economic growth with cycles advanced by 
Goodwin |]2|,|3[] and one of its specific variations, the Desai-Henry-Mosley-Pemberton (DHMP) extension [9], using 
Brazilian income data from 1981 to 2009. The variables used by Goodwin in his model, the workers' share of total 
production u and employment rate v were obtained by describing the individual income distribution by the Gompertz- 
Pareto distribution (GPD) 12711 . formed by the combination of the Gompertz curve, representing the overwhelming 
majority of the population (~99%), with the Pareto power law, representing the tiny richest part (~ 1%) [26]. We 
identified the Gompertzian part of the distribution with the workers and the Paretian component with the class of 
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Unemployment rate [V] (%) 



In (95-[u]) 



Figure 6: Left: [h]/[h] vs. [V]. Right: [v]/[v] vs. [U] (see equationsll4land ll5t . where u was assumed as 95%. In the left plot the full line 
indicates a power-law fit of the form [u]/[u] = A] + B^V] 6 . The fitting parameters yielded A, = -0.011 ± 0.022, B[ = 0.0003 ± 0.0027 and 
S = 1.59 ± 2.94, results which should be compared with equation U2\ . For the right graph the full line indicates a straight line fit using the 
expression [v]/[v] = A2 + B2 ln(95 - [«]). The resulting fitted parameters are A2 = 0.12 ± 0.05 and B2 = -0.054 ± 0.024. This result should be 
compared with equation 113) . 

capitalists and used GPD parameters obtained for each year in the studied time period to analyze the time evolution of 
these variables by means of the Goodwin dynamics. Unemployment data was also obtained from income distribution 
so that all variables come from the same sample since Brazilian unemployment data was collected under different 
methodologies during the time span analyzed here. 

The results were, however, mixed, both qualitatively and quantitatively. The data showed clockwise cycles in the 
u-v phase space in agreement with the model, but those cycles were only largely clockwise and the orbital center was 
not unique, results which brought only partial qualitative agreement of the model with Brazilian data. We obtained 
temporal variations of the variables and their derivatives and carried out straight line fittings to the points formed with 
these quantities, both in the original Goodwin model and its DHMP extension in order to obtain fitting parameters 
which were compared with predictions of both models. In this respect the original model was able to provide a better 
empirical consistency, but the observed parameters were different from what the model predicts in the sense of their 
general behavior, leading to fitted lines whose slopes had opposite behavior than the theory states. A similar situation 
occurred with the DHMP extension, but in this case the uncertainties in the fitted parameters were too large, leading 
to mostly inconclusive results. Although a general predator-prey like behavior was observed, the lack of a single 
orbital center and parameters behaving very differently from what was anticipated bring into question the economic 
hypotheses used by Goodwin in deriving his model. It appears that they may be inapplicable to the economic system 
under study, a conclusion which comes as no surprise in view of the extremely simple specifications of the model, as 
discussed in Sect. 12.41 

Considering these results, in order to provide a viable representation of the real world the Goodwin model must be 
modified. Firstly, as it is obvious from our results, as well as the ones obtained by previous authors, there cannot be a 
single orbital center. We can envisage two possible reasons for such a result: ( i) the "constants" of the model may not 
be constants at all, but time variables; (ii) the right-hand side of equations (O and (0 are too simple and may require 
more terms involving the two variables, which means giving up the linear approximation of equations (f23) and (129V 
In other words, going to a fully nonlinear modeling. 

Secondly, the emphasis so far given by several studies on the economic foundations of the model, which have 
been the main source for its proposed theoretical modifications, should be put aside, at least temporarily, in favor of 
devising differential equations capable of reproducing the observed features, like the moving orbital centers and the 
behavior of the graphs with the temporal variations of u and v. Clearly those economic hypotheses will need to be 
revised as they produce a model which does not agree with the data, but these revisions must be made in the light 
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of empirical results and not solely by theoretical introspection. Possibly new variables representing other economic 
players, like debt and government policy, may have to be introduced in the model, which means that, perhaps, more 
than two coupled differential equations would be necessary to define the economic system. In this respect, as discussed 
by Keen f_8ll . Hudson f4lll and Hudson & Bezemer f4_ll . investment is not profit, being debt-financed when it exceeds 
profit, and government taxation has to be deduced from output to determine profit. 

Thirdly, since the DHMP model fared much more poorly as compared to the original model, Occam's razor 
dictates that these modifications must be focused in the latter rather than the former because the original model is 
simpler. So, developing more complex models without a clear empirical motivation, and in the absence of a clear 
guidance given by real-data observations, goes against Occam's razor. 

The basic motivation behind these proposed modifications comes from the realization that in its present state the 
Goodwin model does not provide much more explanatory power beyond the original qualitative ideas advanced by 
Marx. This is so because it is essentially a mathematical dressing of Marxian ideas by means of a predator-prey set of 
first order differential equations, but which produces solutions that clearly contradict empirical data in many respects 
and provides only general qualitative agreement with real-world observations. Therefore, the real challenge lies in 
devising a model that addresses real-world data and is capable of surviving empirical verification. One must always 
keep in mind that the good scientific practice entails a permanent search of convergence between hypotheses and 
evidences. 
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